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We provide a numerical method to calculate comprehensively the microwave and the laser spectra 
of ultracold bosonic atoms in optical lattices at finite temperatures. Our formulation is built up 
with the sum rules, up to the second order, derived from the general principle of spectroscopy. 
The sum rule approach allows us to discuss the physical origins of a spectral peak shift and also a 
peak broadening. We find that a spectral broadening of superfluid atoms can be determined from 
number fluctuations of atoms, while that of normal-state atoms is mainly attributed to quantum 
fluctuations resulting from hopping of atoms. To calculate spectra at finite temperatures, based 
on the sum rule approach, we provide a two-mode approximation assuming that spectra of the 
superfluid and normal state atoms can be calculated separately. Our method can properly deal 
with multi-peak structures of spectra resulting from thermal fluctuations and also coexisting of 
the superfluid and the normal states. By combining the two-mode approximation with a finite 
temperature Gutzwiller approximation, we calculate spectra at finite temperatures by considering 
realistic systems, and the calculated spectra show nice agreements with those in experiments. 
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I. INTRODUCTION 

Ultracold atoms in an optical lattice allows us to sim¬ 
ulate quantum phase transitions of lattice fermions and 
also bosons [lll3- In fact, the superfluid (SF) to the 
Mott insulator (MI) transition of bosonic atoms has been 
demonstrated by using various measurement techniques 
[Ml- The signature of the phase transitions can be ob¬ 
served in certain thermodynamic quantities il. One 
of examples is to characterize the transition by observing 
the disappearance of a coherent peak structure in the 
number distribution of atoms in the momentum space 
A spectroscopic measurement is another useful 
tool to detect phase transitions fTol - [l5l| . This is because 
much information is included in spectra that reflect the 
dynamical response of many-body systems after excita¬ 
tions caused by a certain external field. Furthermore, 
when the external field is very weak and perturbative, 
the dynamical response can be connected to thermody¬ 
namic quantities of thermal equilibrium states before ex¬ 
citations. In condensed matter physics, such a relation¬ 
ship, e.g., fluctuation-dissipation theorem, has been used 
to discuss quantum many-body phenomena. It is thus re¬ 
quired to deeply discuss such spectroscopic relationships 
specific to cold-atom systems. 

One of pioneering studies on spectroscopic measure¬ 
ment on atoms in a lattice is microwave spectroscopy 
experiments, where the Mott shell structure has been 
observed by spectroscopically distinguishing the differ¬ 
ent number states of atoms [^. Theoretically, the cor¬ 
responding spectra have been studied with an approxi¬ 
mation satisfying the (first order) spectral sum rule di¬ 
ll^, which is derived from the general principle of spec¬ 
troscopy M- The first-order sum rule determines the 
relationship between the spectral peak position and the 
two-body correlation function of atoms fl^ . This is a 


prominent example that connects thermodynamics to dy¬ 
namics in cold atom systems. This calculation assumed 
that the system is at zero temperature, while the real¬ 
istic experiments have been done at low but finite tem¬ 
peratures. In addition, such a first order approximation 
is insufficient to discuss important properties of spectra, 
such as, a standard deviation and a spectral broaden¬ 
ing, which can be connected to fluctuations of atoms in 
thermal equilibrium. On the other hand, the laser spec¬ 
troscopy is now being established [20l - l^ . The laser and 
the microwave spectroscopy are understood as a simi¬ 
lar type spectroscopy based on electromagnet-field ex¬ 
citations. However, the laser spectroscopy cannot be 
straightforwardly described by the formulation of the pre¬ 
vious studies [Ig. This is mainly due to the difference 
in wavelengths of the external fields. A reliable theoret¬ 
ical method for comprehensively analyzing these spec¬ 
troscopy at finite temperatures is now required. 

In this paper, we theoretically discuss a common for¬ 
mulation for the microwave and the laser spectroscopy 
of ultracold bosonic atoms in a three-dimensional opti¬ 
cal lattice. We start with analy zing the sum rules in the 
same way as the previous work [16| , while we extend the 
approximation to the second order. This approach al¬ 
lows us to clarify that number fluctuations of atoms in 
thermal equilibrium can be connected to a broadening 
of spectra. Phenomenological discussions based on the 
sum rule approach allow us to establish a method for 
calculating spectra at finite temperatures. We propose 
a two-mode approximation assuming that the spectra of 
condensed SF atoms and uncondensed normal state (NS) 
atoms are separately dealt with. The multi-peak struc¬ 
tures resulting from thermal fluctuations and also the 
coexisting of the SF and NS atoms can be appropriately 
taken account. Using this approximation combined with 
a finite temperature Gutzwiller approximation we 
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numerically calculate the microwave and the laser spectra 
by considering realistic experimental parameters [111 . |24| . 
We find that our approximations reproduce essential fea¬ 
tures of spectra seen in the microwave experiments [HI, 
and we predict spectra of the realistic laser spectroscopy 
experiments [^ . 

II. THEORY OF SPECTROSCOPY 

This section is devoted to the general theoretical 
framework of spectroscopy. We first explain the model 
Hamiltonian, and then, we show the common formulation 
to describe the microwave and the laser spectroscopy. To 
capture essence of the present spectroscopy, we discuss 
physical properties of spectra in simple model cases. For 
simplicity, we set h = 1 and fcs = 1. 


angular frequency, and a wavevector of the external field, 
respectively, and t and r are time and position. Here Oex 
is an excitation operator defined as follows. For conve¬ 
nience, we define PaOex.a as a second quantization 
of and 

Oex,a = e -|- H.C., (2) 

i 

where Cea,i is the annihilation operator of an excited 
atom in the a-th orbital at the i th site of the posi¬ 
tion Yi. Note that kex is a reduced wavevector in the 
first Brillouin zone defined by kgx = Kgx + G, where G 
represents any reciprocal vectors with e**^''’* = 1. An 
excitation matrix pa is defined by 

p^= (3) 


A. Model Hamiltonian 


Before spectroscopic excitations, thermal equilibrium 
properties of atoms in an optical lattice is well described 
by the following single-band Bose-Hubbard Hamiltonian 


'^g ~ Cg jCg,j + 

<i.i> * 

+ ( 1 ) 

i 

where ^ (cg,i) is the creation (annihilation) operator of 
an unexcited atom at the ith site, and fig^i is the corre¬ 
sponding number operator. Here, contributions of higher 
orbitals can be neglected when we consider the low en¬ 
ergy properties. We note that higher orbitals have a role 
in the spectroscopy as discussed later soon. The Hub¬ 
bard parameters, i.e., the interaction strength f/g g and 
the hopping integral Jg, are evaluated by the ab initio 
calculations based on the second quantization using ex¬ 
perimental parameters: a lattice constant ol and a lat¬ 
tice depth Vo, which are determined from a wavelength 
and an intensity of the lattice laser, respectively, and Og^g 
a scattering length between two unexcited atoms. The 
chemical potential p is determined so as to fix the total 
number of atoms iVtot = J2iAg,i) ^g,i the trap¬ 
ping potential. In the following, for simplicity, we omit 
to explicitly write down p, which can be included in a 
global shift of Vg_i. 


B. Spectroscopy 

In the microwave and the laser spectroscopy, excitation 
processes caused by an external electro- mag netic field are 
generally described by [I^, where F, w, 

and Kex are a non-dimensional normalized amplitude, an 


where WQ,(r — r^) is the a-th Wannier orbital at the ith 
site for the excited atoms, and IFi(r — r^) is that for the 
unexcited atoms (i.e., a = 1). Here, \pa\^ represents the 
probability that the orbital of atoms changes from the 
lowest to the a-th orbital during excitations. The or¬ 
thogonality of the Wannier orbitals assures a condition 
\PoA' = 1- In Ihi® paper, we neglect the probabil¬ 
ity that atoms are excited to the different lattice sites 
(i.e., inter-site excitation), because it is exponentially 
smaller than that of the onsite excitations. Namely, we 
assume that pA = /(irW'*(r — — r^) 

vanishes except for i = j. 

We focus on the weak excitation limit under the con¬ 
dition of |r| <C 1. The excitation spectra can be formally 
given by I{uj) = Ya l/OaP-Ia(w), and 

Ulo) = |(n'|6ex.a|n)|2e-('^'‘-^)/^J(cc-F;;+£;„), 

n' ,n 

( 4 ) 

where \n) is the eigenstate of Hamiltonian T-Lg in Eq. ([1]) 
with energy En, and H(= —Tln^^ is the grand 

potential. Note that the conservation law of the number 
of excited atoms allows us to decompose /(w) into the 
sum of Ia{u!)- Th® excited state |n') is the eigenstate of 
Hamiltonian H = Hg + He + Hge, and E'^ is its energy. 
Here He and Hge are given by 


^ ^ jCea,_7 H” ^ “1“ ^ea,i 


i,Oi 


^ ^ Ug^ea'^ecxA^gA^ 

i,a 

(6) 


where Jea is the hopping integral of excited atoms in the 
a-th orbital, and t/g,ea is the onsite interaction between 
the first orbital unexcited and the a th orbital excited 
atoms. Note that the interaction between two excited 
atoms Uea,ei3 Can be reasonably neglected in the limit of 
weak excitations. Aea represents the energy difference 
between the unexcited atoms in the lowest orbital and 
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the excited atoms in the a th orbital. We can always 
set Aei = 0 by appropriately choosing the origin of the 
spectral frequency. The spectral intensity proportional to 
|rp is determined so as to satisfy the integral condition 
/ I{ui)duj = const., when we compare our analyses with 
the experimental observations. Thus, we can neglect the 
quantitative aspect of T by setting T = 1 without loss of 
generality. 

C. Sum rules 

We discuss the moment expansions of the spectral 
function by following the previous studies [la - ll^ . In 
general, spectra given by Eq. (jT]) should satisfy the fol¬ 
lowing sum rules in terms of the n-th order mo¬ 

ment: 

= J dujuj'^Uuj), 

= ([[[OL.a,^],^r--]Oex.a), ( 7 ) 

where 77], H, • • • ] denotes the n-times commuta¬ 

tor between Oex,a and 'H. These relations indicate that 
certain statistic quantities in thermal equilibrium have 
a relation with some properties of spectral functions re¬ 
flecting a dynamical response of the system. For exam¬ 
ples, by considering up to the second order moments, the 
spectral mean value uja is given by Coa = Ma'^/Ma\ 


and the standard deviation ctq, can be written as cr^ = 
Ma^ /Ma^ — /Ma'')'^. Naively, we can stress that 

Wq and {Tq, determine a spectral peak position and its 
width, respectively, and the relations defined by Eq. © 
can describe the physical origin of the peak shift and also 
broadening caused by the many-body effects. In what fol¬ 
lows, to discuss these important spectral properties, we 
analyze the sum rules up to the second order. 

We can derive the following expressions with respect to 
the corresponding sum rules. The zeroth order is given 
by 

= (OL.aOex.a) = (8) 

i 

The first order moment is written as 
= ([OL,„,77]Oex,o) 

= iUg,ea ~ Ug^g) ^ ^ .i^g.i^g,i^g,i) 
i 

i 

+ ( 9 ) 

i d 

where represents a summation over the adjacent 
sites, and the second order is 


= ([[0|x,a>77],77]6ex.c) 

= {Ug^ea ~ ^^(^g,i^g,i^g,i^g,i^g,i^g,i} + {Ug,ea ~ Ug^g) ^ X ^a.i^a.i^g,i^g,i) 

i i 

H” ^ ^ea) ('^i) ^ ^ ^{Ug^ea 

i i 

+ ‘^{Jg ~ Jea^ ){Ug,eoL — ^g 

i d 

+ E E - ^5,7 + ^eo.){^g/Cg^^+d) 

i d 

+EE {Jg - Jec.e^^- '^){Jg - X (1 ,o){^gd^g,i+d+d') 


i d,d‘' 




( 10 ) 


i d 


The expressions in Eqs. ©-(uni) include the onsite 
multi-body correlation functions = ((Cg_i)^(cg,i)^), 
where Gi,i is equivalent to the averaged number of atoms 
Ui = The hopping Hamiltonian yields the inter¬ 

site correlations such as Qij = (c^ jCgj), where Qi^i = Ui, 


and = (Cg_iCg_icVcgj). In the second (or higher) 

order of moment, the hopping Hamiltonian yields also 
the onsite one-body correlation Gi_i(= Ui) [see the last 
term in Eq. (nni)]. Note that this type of terms is caused 
by a round-trip hopping process. These terms are pro¬ 
portional to Jg — where describes the 
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momentum transfer from the external field to atoms. The 
momentum transfer can be regarded as a back action of 
the measurements, which plays an important role in spec¬ 
tra. 


D. Physical meaning of the spectral deviations and 
spectral mean value 

It is convenient to discuss the physical meaning of 
spectra, which can be figured out from the sum-rule 
approach. The spectral mean value Ua can be writ¬ 
ten as the sum of 5Uon and 5Ja^ which are 

spectral energy shifts caused by the effects of interac¬ 
tion, trapping potential, band gap, and hopping, re¬ 
spectively. As discussed in the previous study 0, 
5Ua = {Ug,ea “ Ug,g) ^ 2,1 /Ntot is a collisional en¬ 

ergy shift. Two terms, SVa = J^i(Vea,i - Vg^i)ni/Ntot 
and 5/\a = arc related to the statis¬ 

tical average of the number of atoms. Here, (5 Aq, re¬ 
duces to a constant Ago, having no connection to any 
thermodynamic quantities, while 5Va describes the ef¬ 
fects of the inhomogeneity of the system. The hopping 
energy shift ~ Ntot 

is related to the inter-site correlation. The Bloch band 
picture makes physical meanings of this term clear. We 
define the kinetic energy shift 5Ka = 5 Ja SAa, which 
can be rewritten as X]k(^ea,k+ke„ -«?s,k)%,k/Atot, where 
7 is the momentum space distri¬ 

bution, and Eg.k and Ega.k are the dispersions of the unex¬ 
cited atoms in the lowest orbital and of the excited atoms 
in the a-orbital, respectively. Because of the momentum 
conservation law, a wavevector kex is transfered from an 
external field to atoms. It turns out that 5Ka becomes 
large for a large kgx- Note that the discrete translational 
symmetry imposes kgx to be a reduced wavevector in the 
first Brillouin zone. 

We next discuss the standard deviation Ua related 
to the second order moment as cr^ = — 

{Ma ^Roughly speaking, the standard devi¬ 
ation characterizes broadening of spectra. For clear 
vision, we now focus on a uniform system by setting 
Vea,i — Vg^i = 0. Considering the physical origin of the de¬ 
viation, we can rewrite a a as = ctu ^ + cr^ „ -I- , 

' UL c/,a: rv,u: y/UK^Ot' 

where the deviation induced by the correlations au,a is 
defined as = (Ug^ea - C^g.g)^Ei(G'3,i + G 2 ,i)/Ntot - 
G' 2 ,i/Atot)^], and that caused by kinetic terms ax,a 
is defined as = I]k(^g,k - eea,k+ke,)^ng,k/Atot - 
Ek(^eQ,k-i-ke,, - eg,k)ng,k/Atot]^, and that originating 
from the cross terms a^ is given by ^ = 

^Ug^ea - C^g,g)Ed('^g “ Ei “ 

(Si G’2,i/Atot)(Si Si,i+d/Aftot)]- 

In the following, we consider two simplified model cases 
and calculate a mean value Qa and a deviation (Tq,, so as 
to discuss what are the physical origins of spectral peak 
shifts and broadening. Here we focus on uniform systems 
at zero temperature for simplicity, and inhomogeneous 


systems at finite temperatures will be discussed in Sec. 
imi We first consider SF states at zero temperature and 
use the following simple mean-field approximation. We 
assume that (Cg,i) is finite and is given by a classical com¬ 
plex number c^, which leads to ^ (CgiCgj) ^ 

c*Cg, and rig^k ~ i^k.koAsp, where A’sf is the number of 
the condensed SF atoms, and ko is a wavevector at the 
bottom of a dispersion eg ^ (usually kg = 0 for a positive 
Jg). We here also assume that almost all atoms are in the 
condensed state A’sf ^ Atot. The spectral mean value is 
now given by = (I/„.g„ -Un.n)ni + (egg,.kn+k., - £».kn). 
It is consistent with the previous study [1^, while the 
additional kinetic energy shift is found. The deviation 
now reduces to cr^ = {Ug^ea - Ug^g)‘^ni{= crfj^^). In¬ 
terestingly, we find that the deviation has a contribu¬ 
tion of an interaction term only. Here because of 

the subtraction in Ma'^/Ma'^ — [M^'^/Ma'^Y , the de¬ 
viations aK,a and a ^ are canceled out under the 
pure-condensation condition rig k ^ i^k.ko Atot- The num¬ 
ber of the coherent SF atoms is indefinite, and number 
fluctuations An^ of such a coherent state are written as 
y/ni, where {Am)"^ = (hg J — nf. Thus, we can conclude 
that the spectral deviation of coherent SF states is con¬ 
nected to number fluctuations in thermal equilibrium as 

— \Ug^ea ^gjglAll^. 

We next consider MI states with m atoms in each 
site at zero temperature. For MI states, it is reason¬ 
able to set inter-site correlations Qij for i ^ j to be 
zero. The spectral mean value is now given by Wq, = 
{Ug,ea — Ug^g)(rn — 1) -|- Ago,, which is equivalent to the 
previous study [l^ . The spectral deviation reduces to 

= J2d l-^g “ Jeae*’^“'‘^p(= In the same man¬ 

ner as the above, the correlation-induced deviation au,a 
cancels out as ufj^ = m{m—l){m—2)/m+m{m—l)/m— 
[m(rn — = 0. The cross-term-induced deviation 

'^VuK a because of negligible inter-site corre¬ 

lations. For MI atoms at zero temperature, the number 
of atoms is definite, and there are no number fluctua¬ 
tions Arii = 0. On the other hand, the phase and the 
momentum is indefinite, and thus a spectral broadening 
of MI atoms is caused by kinetic fluctuations described 
by o'K,a- By executing the k-summation in <JK,a with a 
constant momentum distribution ng^k(= m) of uniform 
MI states, we obtain ^ | Jg — Jeo,e®*'“''^p. This 

deviation can be connected to quantum fluctuations re¬ 
sulting from the round-trip hopping process given in the 
last term in Eq. (ITO)) . 

As demonstrated by the previous experiments [TTj | , the 
present spectroscopy has an ability to distinguish the dif¬ 
ferent number states, when |17g^ea “C^.gl large enough. 
This feature of spectra can be explained by the first or¬ 
der sum rule, Wa = (Un.pn — Ug^g)(m — 1), as discussed 
in the previous study [l6|. The above second-order sum 
rule approach further indicates that number fluctuations 
of atoms in thermal equilibrium An^ play an important 
role in such number-resolving spectroscopy. In fact, a 
spectral deviation ctq, of coherent SF atoms can be deter- 
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mined from number fluctuations: ctq, = \Ug^ea — Ug^g\Ani. 
Even for MI atoms with An^ = 0, a spectral deviation 
is finite owing to kinetic fluctuations, which is attributed 
to the indefinite phase and momentum in a reflection of 
the definite number and position of MI atoms. This de¬ 
viation of the number definite states is not related to 
any thermodynamic quantities, and thus this constant 
ctq, of ~ is the intrinsic lower limit 

of a spectral linewidth (see Sec. IIII Bl) . On the other 
hand, the deviations of general states with both phase 
and number fluctuations are given by the summation, 
’^u,a + '^K,a + '^^/uK a' Spectral de¬ 

viations of a specific ground state (SF and MI) with a 
definite quantity (phase and number) are characterized 
by fluctuations resulting from the conjugate indefinite 
quantity (number and phase, respectively). It should be 
noted that the second-order sum rule makes clear the 
fact that the spectral measurements are governed by the 
uncertainty principle. 

These physical properties of spectra mentioned above 
are figured out from the general principle of weak exci¬ 
tation spectroscopy with given Hamiltonian 'H and an 
operator Oex,a- For example, an excitation operator 
Oex,a characterizes the intrinsic spectral broadening of 
\ Jg — where a momentum transfer of 

kcx, which is a back action of the measurements, deter¬ 
mines a quantitative aspect of a spectral width. Bose- 
Hubbard Hamiltonian Tig includes kinetic and interac¬ 
tion terms, and the competition between these two con¬ 
jugate terms is the origin of the SF-MI transitions. The 
second-order sum rule approach clarifies that spectral 
deviations reflect the completely different properties of 
these two conjugate states, SF and MI. We can thus 
conclude that the present spectroscopy will be a sensi¬ 
tive tool for detecting the SF-MI transitions. Note that 
the first order sum rule approach is insufficient to clarify 
these important features of spectroscopy. However, the 
above simplified discussions cannot be straightforward 
applied to the finite temperature spectra. In the next 
section, we thus propose a two-mode approximation to 
numerically calculate spectra that satisfy the sum rules. 


III. METHODS 


In this section, we provide a numerical method for 
calculating spectra at finite temperatures. We first ex¬ 
plain the finite temperature Gutzwiller approximation 
1^ . which allows us to efficiently obtain the thermo¬ 
dynamic quantities in Eqs. (l5])- (fTUl) . We next provide a 
two-mode approximation to numerically calculate finite 
temperature spectra in inhomogeneous systems. At the 
end of this section, we compare our method with the pre¬ 
vious formulations BUS- 


A. Finite temperature Gutzwiller approximation 


The Gutzwiller approximation allows us to efficiently 
analyze the thermal equilibrium properties described by 
the Bose-Hubbard Hamiltonian in Eq. ([T]). This is a 
mean-field approximation considering up to the first- 
order collection in terms of Jg and well describes the 
SF-MI transitions in high dimensional systems. Here, 
Jig is then approximated by a set of the effective local 
Hamiltonian Hioc = with 




efF* ; 
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1 ); 
( 11 ) 

where J®® is determined from a self-consistent condition 


= —J2d’^gi^g,i+d)- Using exact diagonalization, we 
can numerically calculate statistical quantities such as 
Ci = (Cg^i) at finite temperatures [^ . 

As discussed in Sec. IHDI a finite Ci effectively describes 
the Bose-Einstein condensates (BEG) within the mean- 
field approximation. Here, the number of condensed SF 
atoms in each site can be defined by usf,z = \ci\^■ Both 
thermal fluctuations and interactions cause coexisting of 
condensed SF and uncondensed NS such as MI and nor¬ 
mal fluid (NF). The annihilation operator of NS atoms 
at the fth site is effectively given by CNS.i = — Ci, 

where NS atoms satisfy always a condition (cns.i) = 0- 
The number of NS atoms riNS.i = (c^g ffiNS,*) can be 
written as nNS,i = ni — nsF,i- This leads to the fol¬ 
lowing reasonable condition: Wot = Nsf + A^nSj where 
-^SF(NS) = X)WSF(NS),i- The total number of atoms is 
given by the sum of the total number of SF and NS atoms. 

It is useful to briefly explain how to calculate the ther¬ 
mal quantities in the moments Eqs. ([9]) and (HOD . The on¬ 
site correlation functions Gn,i can be calculated straight¬ 
forwardly by diagonalizing the effective local Hamilto¬ 
nian 'Hioc,i- On the basis of this local approximation, the 
inter-site correlation is described by Qij ^ (c^ i)(cgj)(= 
c*Cj) for i j. The higher-order inter-site correlations 
G 2 ,i,i+d also reduce to (Cg ^CgyCg J(cgj). These expres¬ 
sions mean that the inter-site correlations of NS atoms 
(cjjs ffiNS,i) for i ^ j are approximately set to be zero, 
and thus, two kinds of NS states, MI and NF, are dealt 
with approximately in the same way. We should note 
that the MI states appearing at lower temperatures can 
be characterized by focusing on the creation of the Mott 
shell structures and also the suppressed entropy per site 
[ 2 ^. We can thus effectively calculate that thermal fluc¬ 
tuations cause the MI-NF crossover within this local ap¬ 
proximation. 


B. Two-mode approximation 

Next, we provide a two-mode approximation that helps 
us to calculate spectra at finite temperatures. We assume 
that Ia{oj) in Eq. dH) can be decomposed into two com- 
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ponents resulting from the contributions of SF and NS 
atoms: 

= (12) 

Two types of uncondensed NS (MI and NF) states ap¬ 
pear at finite temperatures. As mentioned in Sec. IIII Al 
within the Gutzwiller approximation, NF states are ap¬ 
proximately dealt with in the same way as MI states 
based on the local Hamiltonian picture. As discussed 
in Sec. Ill PI at zero temperature, spectra of coherent SF 
atoms show completely different properties by comparing 
with those of MI atoms. Note that the special charac¬ 
teristics of spectra of SF atoms result from the phase 
coherence caused by BEG. We thus deal with spectra of 
SF atoms in different way to two types of NS atoms. 

On the basis of the two-mode approximation, we re¬ 
consider the sum rules for the spectral moments: 

= J dujuj^l^^{uj)+ j dww^lf{u;). 

The sum rules up to the second order (i.e., up to n = 2) 
provide the following relations: 

= A^sf + A’ns, (13) 

mW = iVj,s^NS ^ -SF^ (^4) 

= [(uf)" + + (caf )2]iVsF, 

(15) 

where and are the spectral mean value, and uns 
and (TsF are the spectral standard deviation for the NS 
and SF spectra, respectively. The zeroth order sum rule 
in Eq. (fT^ simply offers the condition associated with the 
total number of atoms, which is always satisfied within 
the Gutzwiller treatment as mentioned in Sec. Hira On 
the other hand, the first and second order sum rules re¬ 
quire the balance conditions between /^®(a;) and /®^(a;), 
and these conditions allow us to properly calculate spec¬ 
tra. 


1. Spectra of uncondensed normal state atoms 


The spectral weight Wa,i,m and the peak position Pa,i,m 
are given by 

Wa,i,m = (17) 

Pa,i,m — l)(G^^ea ^g,g) ^ea “b ^ea,i ^gd^ 

(18) 

where is the Boltzmann factor of the 

number state \m)i, and Hi(= —Tln^^ is 

the grand potential in the i th site. Note that, even 
for the uniform systems, spectra at finite temperatures 
have multi-peak structures depending on the thermal 
distributions of the number states \m}i described by 

We now discuss that the zeroth and the first or¬ 
der sum rules are always satisfied in the above expres¬ 
sion in Eq. (fTHll when Qij <C u-i. By using Eqs. (HID 
and (IT51) . we obtain Ma'^ = m 

and - l)(Clg.ea - Ug,g) + 

Aea + 14a,z “ Using relations 

and G 2 ,i = ~ 

l)g-/ 3 (Sm,i-Oi)^ we find that Ma'^ = A 4 is, and re¬ 
duces to ^j[(Ug^ea Ug^g^G 2 ^i (A eoL “ 1 “ 

These facts suggest that in Eq. (HU) reproduces 

the zeroth and the first order moments in Eqs. ([ 8 ]) and 
(|n]) when we can neglect the term proportional to the 
inter-site correlations [the last term in Eq. (O]. 

In contrast, the second order sum rule is not straight¬ 
forward. When Qij <C Ui, almost all terms are repro- 

f 2) 

duced in the same way as the above. Namely, Ma = 
Si m 'uta,i,mP^ i m 1 ® equivalent to the first four terms 
in Eq. (HHl). However, we cannot reproduce one of 
the terms in Eq. (HOI) , which is the round-trip hopping 
term given by Si Sd I'^g ~ Jeae*’^“'‘^pni. This means 
that, even though the inter-site correlations are negli¬ 
gible, quantum fluctuations resulting from the round- 
trip hopping broaden the spectral width of each peak 
in Eq. (fTBI) . Namely, the sum rule requires that the delta 
function S{uj) in Eq. (fTHl) should be replaced with a cer¬ 
tain function with a finite spectral width. We here use a 
Gaussian function, and /„®(a;) is now given by 


In what follows, we discuss the properties of /^®(w) 
and /®^(a;) at finite temperatures, separately. Here we 
begin with J^®(a;) by assuming that Nqf = 0 , and ac¬ 
cordingly J®^(a;) vanishes. We also assume that the inter¬ 
site correlations are negligible Qij for i ^ j hy comparing 
to the on-site correlations n^. The density matrix of such 
a localized state is given by 

at finite temperatures, where Em,i = Ug^gmirn — l)/2 — 
Vg^iTTi is the energy of the local number state \m)i. The 
spectra can be obtained in a form of the exact represen¬ 
tation: 

(^) ~ ^ Padm^' (^^) 

i,m 


exp {-{uj - / {2yl)) 

7a 

(19) 

where 7 a is the spectral width defined by 7 ^ = | Jg — 

Jeae*’^“'‘^p. We comment that a Lorentzian func¬ 
tion is not suitable for the substituting function, be¬ 
cause the second order moment does not converge: 
-|- 7 ^) —>■ 00 . 

The extended representation in Eq. (IT^ with Eqs. (flTl) 
and m properly satisfies the sum rules up to the 
second order when Qij rii. In the same way 
as the above, we can straightforwardly confirm that 


/a"(^) = 


E 


w. 
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the zeroth and the first order sum rules are satis- 
fied. The second order Ma is extended as fol¬ 
lows: J |{■^aV^)duJ = 

E*,m The first term 

'^a,i,niPa^i^rn equivalent to the second order 
moment obtained from the original representation in 
Eq. (dni). The additional term 7 ^ J2i,m Wa,i,m{= 7 a-^Ns) 
properly describes the last term in Eq. (HOD . 

We here estimate the magnitude of 7 a that character¬ 
izes an intrinsic spectral broadening caused by hopping- 
induced quantum fluctuations. For simplicity, we con¬ 
sider a = 1 and set Jg ^ Jei, which leads to 7 ^ = 
2 z [1 — ^ j cos(kex ■d)/z\, where z is the number of the 
neighboring lattice sites (z = 6 in the cubic lattice). For 
kcx ~ 0, 7 i reduces to zero. For kex ~ ( 7 r, 7 r, 7 r), 71 takes 
a maximum 2 Y^|dei| = VFi/i/z, where Wa = 2 z|Jea| 
is a bandwidth of the a th orbital. We next consider 
higher orbitals a 7 ^ 1 by assuming \Jg\ <C | Jea\ , and then 
we obtain 7 a = \/z\Jea\ = Wal‘2-\/z- Simply put, the 
kinetic spectral broadening is proportional to the band¬ 
width, 7 a oc Wen where the wavevector conservation law 
determines the proportionality coefficient ranging from 0 
to \fz depending on kex- 

Before closing the discussions on we consider 

the validity of the condition Qi^j <C n^. For = 0, this 
condition is exactly satisfied: Qi^j = 0 for * 7 ^ j. For a 
finite but small Jg(<C C7g_g), where MI states will appear 
at low temperatures, the effects of interactions strongly 
suppress the inter-site correlations. Large potential dif¬ 
ferences strongly suppress the inter-site correlations (e.g., 
\Vg^i — Vg^j \ » Jg), und thormal fluctuations also decrease 
Qi^j. We thus expect that the NS atoms in the realistic 
systems with interactions and trapping potential at fi¬ 
nite temperatures will satisfy well the condition of small 
rii). This condition is equivalent to flattened mo¬ 
mentum distributions rig.k, which can be confirmed in ex¬ 
periments by using the time-of-flight measurements with 
the projection onto the first Brillouin zone 


2. Spectra of the superfluid atoms 


broadening caused by number fluctuations, which allows 
us to reasonably use this simple assumption. 

To compare /^^(w) with here we mention 

again the properties of at finite temperatures. As 

shown in Eq. (fTO)) . has many peaks at positions 

Pct,i,m(oc Ug^ea — Ug,g) with a Spectral width of7a(oc Wa). 
Note that peak positions and a width are usually deter¬ 
mined from the different energy scales. Thermal fluctu¬ 
ations change relative spectral weights and also increase 
the number of spectral peaks. This multi-peak structure 
is an essential feature of /a®(w) at finite temperatures, 
which can be properly dealt with in Eq. (I19D . In contrast, 
as discussed in Sec. imi for SF atoms, a spectral mean 
value is proportional to the average number of atoms 
oc (Ug^ea — Ug^g)n, and a spectral deviation is given 
by Wg,ea — Ug^g\\/fL. This fact suggests that 

and tr®^ will be usually comparable. A main role of ther¬ 
mal fluctuations is a decrease in the number of condensed 
atoms A^sf- We thus expect that the essence of /|^(w) 
at finite temperatures can be captured by a single peak 
broadened by large number fluctuations. The effects of 
decrease in A'sf on the spectra are considered via calcula¬ 
tions on the sum rules in Eq. (TT^ . and also Eqs. m and 
(TTf)]) . which properly describe decreases in a peak height, 
a peak shift and broadening, respectively. 


3. Spectra for coexisting region 

We next explain the formulation for the middle region, 
A^sf 0 and A^ns 0, where NS and SF atoms coexist. 
Here, /®^(w) and Ia^{uj) are separately calculated, and 
/a®(w) is the first. We assume that NS atoms are affected 
by a mean-field potential resulting from interactions with 
SF atoms. Thus, we here consider the following effective 
local Hamiltonian of NS atoms excluding SF atoms: 

"HnS = bsF,i)^NS.i 

i 

+ Ug,g flT^S,iin-NS,i " l)/2, (21) 


Next, we consider the opposite limit, Nsf ^ -^nSi 
where we neglect Taking account of the phys¬ 

ical properties of BEG, we assume that /a^(w) can be 
described by the following single peak structure: 


= -^SF 


exp(-(a;-:7®nV2( a. 

((T®^V^) 


SF',2 


)^) 


( 20 ) 


where a spectral peak position w®^ and a spectral width 
cr®^ are determined from the sum rules in Eqs. (I14D and 
m, respectively. For such a single peak structure, the 
deviation coincides with the spectral width. This single 
peak assumption may be oversimplification. It should 
be noted that we carefully take account of the spectral 


where the potential VsF,i describes effectively mean-field 
interactions between SF and NS atoms given by Vspy = 
2Ug^gnsF,i + 5Pi- To impose a self-consistent condition 
= {ng,i) - |(cgy)P, we further define a chem¬ 
ical potential shift Spi, where (• • • )g^^ is the statistical 
average at thermal equilibrium definei§ by the Hamilto¬ 
nian ?7 nS: while {■■■) is that defined by the localized 
Hubbard Hamiltonian in Eq. (fTTI) . We note that Spi ~ 0 
for nsF.i ^ ?T-NS,i or nsF.i ^ n-Fis.i, because the mean- 
field treatment is appropriate for these dilute regions. 

We summarize a procedure for calculating the full 
spectra I{uj) = /®^(a;) -I- /^®(w), where /®^(w) = 
Ea and ^ |pap/^®(^). 

1. We first calculate thermal equilibrium states of 
Hubbard Hamiltonian based on the finite temper- 





ature Gutzwiller approximation. We use exact di- 
agonalization to solve the localized Hamiltonian in 
Eq. (fTTll at finite temperatures. We obtain the mo¬ 
ments from Eqs. ©-(Uni), and other statistical 
quantities such as nsF,i and riNS.i- 

2. Next, we calculate in Eq. (IT^ by exactly 

diagonalizing the effective Hamiltonian in Eq. (I^Tl) . 
After that, we can directly calculate and 
from For the self-consistent condition men¬ 

tioned above, we need and Ci, which should be 
obtained in the previous process 1. 

3. Finally, we determine and by using 

Eqs. (fHl) and (fT^ . where we use M^\ (D^®, and 
cr^® obtained in the previous processes 1 and 2. 
Then, we can calculate the full spectra /(w). 

In this way, based on the sum rule approach, completely 
different features of /^®(u;) and /®^(a;) are properly dealt 
with, and multi-peak structures resulting from finite tem¬ 
perature effects and coexisting of SF and NS atoms will 
be taken account of precisely. 

C. Role of the sum rules in the two-mode 
approximation 

We finally discuss what a role the sum rules play in 
the present calculation procedure. The sum-rule ap¬ 
proach combined with the two-mode approximation pro¬ 
vides us with the reasonable relationship between 
and As discussed below, w®^ and cr®^ can be de¬ 

termined reasonably within the level of the mean-field 
treatment. 

We here consider the first order moment Ma\= 
A^NS^NS.a + NsFidsF.a) and focus On the effects of inter¬ 
actions [the first term in Eq. dH])]: {Ug^ea — Ug^g) G 2 ,i, 
and other terms are neglected for clarity. Using a defi¬ 
nition of CMS,i = Cg,i — Ci, we can rewrite G 2 ,i as -I- 

4nNS.i?T-SF.i +nlp^i, where G^f = (cNs.z^NS.iCNS.iCNS.i)- 
Three terms, G^f, nNS.i''^SF,z, and n-gp^, represent the 
two-body correlations between two NS atoms, between 
NS and SF atoms, and between two SF atoms, respec¬ 
tively. As discussed in Sec. lHIBTl the first order moment 
of /^®(w) is easily obtained: 

= {Ug^ea - y^(G^® + 2nNS,i»^SF.i)i 

i 

where the latter term results from the mean field po¬ 
tential VsF,i in Eq. (1^ . Consequently, the sum rule in 
Eq. (1141) allows us to determine the correlation term in 
the first moment of /®^(a;): 

A^SFW®^ = (t/g,ea — Ug^g) y^(n|p ^ -I- 2nsF.inNS,i)- 

i 

The SF-NS correlation 4nNS,insF,z is shared equally be¬ 
tween /^®(a;) and /®^(a;). We now again consider the 


uniform system to compare this expression with those in 
Sec. KDl The collisonal energy shift for the SF atoms 
is now given by w®^ = (Ug^ea - Ug^g){nsF,i + 2nNS,i)- 
The first term (Ug^ea — Ug^g)nsF,i is equivalent to that 
of the pure SF atoms as discussed in Sec. Ill PI and also 
in the previous study [l^ . while the second term is the 
additional contribution originating from the mean-field 
NS-SF interactions. 

In the same way as the above, we can obtain 
the correlation terms in the second order moment of 
/^®(u;) and = (Ug,ea - 

+ 4nsF.i) + 4nNS,znip,J, and 

NSFlia^? + = (C/g,e« - C/g,g)2Ej^iF. A 

^SF,i(l + SuNS.i) + nsF.i(4n^g_j -f 5G^®)]. For the uni¬ 
form system, (cr®^)^ is written as (t/g,ea — Ug^g)^[nsF,i + 
nsF,ims,i + G^® -H 4(AnNS.i)^], where (Att-ns,*)^ = 
(^NS i) ~ ''^NS i- Note that AnNS,i ~ 0 at low tempera¬ 
tures. The first term corresponds to number fluctuations 
of the SF atoms, and the other terms suggest that the 
NS-SF correlations enhance number fluctuations of the 
SF atoms and further broaden /®^(w). 

D. Comparison with the previous studies 

Here, to discuss differeuce betweeu the preseut aud the 
previous treatmeut, we summarize the previous formu¬ 
lation [I^, which has been successfully applied to the 
microwave spectroscopy experiments [111] : 

-f(w) = ^ UiS (w - (C/g,ei - Ug^g)G 2 ,i/ni). (22) 

i 

This approximation satisfies the zeroth order sum rule 
Ma'^ = Wot and partly satisfies the first order, Ma^ = 
J 2 i^ 2 ,iiUg^ea — Ug^g), while the second order sum rule 
is not satisfied at all. Note that, without a considera¬ 
tion of the deviation CTq,, the spectral broadening is ef¬ 
fectively described by the site-dependent G 2 ,i and rii re¬ 
sulting from the inhomogeneity of the system. When 
fcex = 0, j Jg - = 0, Vea,i “ Vg,i = 0, T = 0 and 

the coexisting of NS and SF states is neglected (either 
-Nns = 0 or Asf = 0), the spectral position w®^ obtained 
from Eq. © or Pi^a,m in Eq. (fTSl) and the correspond¬ 
ing term (t/g,ei — Ug^g)G 2 ,i/ni in Eq. (1^^ are equivalent 
with each other. Thus, within the first order approx¬ 
imation, our method is consistent with Eq. (1221) in the 
following two limits; pure SF or MI phases at zero tem¬ 
perature. We note that, at least, the conditions fcex = 0, 
|dg —= 0, and I4a,i —= 0 are well satisfied 
in the microwave spectroscopy (see Sec. II V A")) . 

IV. NUMERICAL SIMULATIONS 

In this section, we show numerical results calculated 
by considering realistic parameters of the following two 
experiments; the microwave-spectroscopy of ®^Rb atoms 
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and the laser spectroscopy of atoms. We first ex¬ 

plain the parameters and then discuss obtained results. 


A. Parameters 


We first point out intrinsic differences of the microwave 
and the laser spectroscopy, and provide parameters used 
in the calculations. A length of microwave is much longer 
than a lattice constant aL(= Al/ 2) and a lattice laser 
wavelength Az, (e.g., of 1064 nm), and as a result, a 
wavevector fcex can be set zero. Therefore, the parameter 
region of the microwave spectroscopy corresponds to the 
perfect Lamb-Dicke regime, where excitation matrices pa 
in Eq. (|31) are given by pi = 1 and Pa/i = 0. Here the ki¬ 
netic energy shift 5Ka and deviations (JK,a and a^ 
are also negligible. In contrast, for the laser spectroscopy, 
a wavelength of the excitation laser Aex (e.g., of 507 nm) 
is comparable to Az (e.g., of 532 nm). Namely, fcex is the 
same order as a lattice wavevector 27r/az. Thus, orbital¬ 
changing excitations and kinetic contributions will be im¬ 
portant in the laser spectroscopy. 

The microwave spectroscopy of ®^Rb atoms uses an 
excitation between different hyperfine states [llj |. while 
the laser spectroscopy of ^^^Yb atoms uses an excitation 
between different electron configurations, ^Sq and ^P 2 
states Rb atoms are trapped by the combina¬ 

tion of optical and magnetic potential, and Yb atoms are 
trapped with optical potential. For the microwave spec¬ 
troscopy, we use a harmonic trapping potential Vgp oc 
Cxxf + Cyuf -I- Czzf, where the curvatures Cx, Cy, and 
Cz are determined from the experimental parameters so 
as to reproduce the bottom of trapping potential till, 
while for the laser spectroscopy, we use an anharmonic 
potential by carefully considering the laser configurations 
in experiments M- Note that the trapping potential of 
two hyperfine states was set to be the nearly same m, 
while ^Sq and ^P 2 states of Yb atoms are trapped in the 
different potential due to the greatly different polarizabil¬ 
ity. Thus, Vea,i — Vg^i Can be set zero for the microwave 
spectroscopy, while it is finite for the laser spectroscopy. 

Differences in scattering lengths ag^e — o-g,g are -0.13 nm 
and -30 nm for the microwave and laser spectroscopy, re¬ 
spectively. Both of them are negative, so that differences 
in the interaction strengths Ug^ea — Ug^g are also negative. 
In addition to the spectral broadening caused by quan¬ 
tum fluctuations (see Sec. IllDl) . we consider linewidths of 
the excitation laser of about 1 kHz and of the microwave 
of about 5 Hz, including the Fourier width of the excita¬ 
tion pulse. 


(a) microwave spectroscopy (b) laser spectroscopy 



-160 -80 0 80 160 -40 -20 0 20 40 


m[Hz] 


m[kHz] 


FIG. 1: (Color Online) Spectra I{uj) for the microwave and 
laser spectroscopy in a deep lattice calculated by the following 
parameters: (a) Vo = 35Er, Ntot = 10®, and T = 100 nK, and 
(b) Vo = 15i?r, Ntot = 2.2 X 10'*, and T = 100 nK. There are 
no SF atoms because of the strong interactions; Nsf 0 and 
/(co) ~ 1^3(0.). 


rameters in the microwave spectroscopy are Vq = ibEx 
and Ntot = 10®, and those in the laser spectroscopy are 
Vb = IbEx and Ntot = 2.2 x 10*, where Er is the recoil 
energy. A spectral peak appearing at w = 0 always means 
that the m = 1 number state (|m = 1)) is excited without 
orbital changing, because the origin of spectra is renor¬ 
malized by setting Agi = 0. Since Ug^ea — Ug^g is negative 
for both spectroscopy, peaks of \m) with m > 2 appear 
orderly in the region of w < 0. The orbital-changing 
excitation requires a large positive bandgap energy Ae 2 . 
Thus, spectra in the laser spectroscopy show some peaks 
in w > 0, which have the similar characteristics to those 
in a; < 0 but have the small intensities because of a small 
excitation probability |p 2 p(~ 0.1|pip). 

Next, we discuss a spectral peak width. For a deep lat¬ 
tice, spectral broadening is mainly attributed to the ef¬ 
fects of inhomogeneity. Equation ([9]) shows that potential 
energy difference (14a,i — Vg^i)ni/Ntot causes just a peak 
shift. However, due to the inhomogeneity, a variation 
of (14a,i — Vg^i)ni/Ntot for different i effectively induces 
broadening of spectra. For the laser spectroscopy, an en¬ 
ergy scale of this variation is estimated to be about 1 kHz, 
which is consistent with the obtained spectral features in 
Fig. [T] In the present parameter region, broadening ef¬ 
fects resulting from the hopping terms are negligible by 
comparing with this inhomogeneous broadening. On the 
other hand, since |14a,i — NgA ~ 0 for the present mi¬ 
crowave spectroscopy, a width of each peak nearly equals 
to a linewidth of the microwave pulse. 


C. Lattice depth dependence of spectra in the 
microwave spectroscopy 


B. Comparison between two spectroscopy for the 
deep lattice 

We first discuss spectra in a deep lattice and compare 
two kinds of spectroscopy. Figure [U shows the spectra 
calculated at a temperature T of 100 nK. The other pa- 


Next, by focusing on the microwave spectroscopy, we 
discuss how spectra change as the lattice depth varies. 
Figure [2] shows the spectra I{uj) calculated with Vq = 
5Ex, lOEx, and 25Ex for Wot = 10® at T = 25 nK. For 
Vo = 25Ex, the number of SF atoms Wf is zero, and a 
discrete peak structure resulting from appears. In 
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FIG. 2: (Color Online) Spectra I{io) of the microwave spec¬ 
troscopy for Ntot = 10® at T = 25 nK for Vo = 25Er, lOEr 
and 5Er. For Vo = WEr, the contributions of SF and NS 
atoms, and are also shown by thick and thin 

lines, respectively, where I (to) = For 5Er, 

A^ns ~ 0 and I{uj) ~ 7®^(cu), while for 25Er, Nsf 0 and 
/(co) ^ /NS(cc). 


contrast, for Vq = 5Er, almost all atoms are in condensed 
states A^sf ~ A^tot and A^ns ~ 0. From discussions in Sec. 
Ill C[ we can naively expect that the spectra of SF atoms 
I^^{uj) are centered at around {Ug^ea — Ug^g)n, and a peak 
width is determined from \Ug^ea — Ug^g\\/E, where n is the 
averaged number of atoms. From these features, we can 
estimate fi ^ 2. For a middle region Vq = 10Er, we 
find an asymmetric broadened peak structure. We note 
that this characteristic asymmetric structure can be at¬ 
tributed to the coexisting of NS and SF atoms. Figure 
m also shows spectra of SF atoms and those of 

NS atoms The spectra of SF atoms are 

centered at around 2.5{Ug^ea — Ug^g). On the other hand, 
has a large intensity at around w ~ 0 correspond¬ 
ing to the \m = 1) excitation. Here we find no discrete 
peak structures, because \ Ug^ea — Ug^g\ is smaller than the 
microwave linewidth. The sum of two spectra yields the 
characteristic asymmetric spectra. 

The obtained spectra for all three parameters in Fig. 
[2] capture essential features of those observed in experi¬ 
ments [HI: The experimental spectra show the symmet¬ 
ric single peak structure for Vq = 5Er, and the asymmet¬ 
ric broadened peak structure, or an overlapped double¬ 
peak structure, for Vq = lOEr, and the discrete peak 
structure for Vq = 25Er. Our assumptions, e.g., the 
single peak assumption for the SF atoms, and the two¬ 
mode approximation, appropriately reproduce character¬ 
istic structures of the spectra seen in experiments. 


D. Temperature dependence of spectra in the 
microwave spectroscopy 

Next, we discuss the temperature dependence of spec¬ 
tra by focusing on the microwave spectroscopy again. In 
Fig.O we first show the spectra for Vq = 35Er at T = 100 



7=100 nK ^ 
50 nK - 
10 nK ■ 


03 [Hz] 


FIG. 3: (Color Online) Spectra /(w) for the microwave spec¬ 
troscopy for Fo = 35Er and Ntot = 10® at T = 100 nK, 50 
nK, and 10 nK. 



FIG. 4: (Color Online) In-situ distribution of the number 
states \m) for Vo = 35Er and Ntot = 10® at T= 100 nK (a) 
and 10 nK (b), and cross sections along the (white) line. 


nK, 50 nK, and 10 nK. From Fig. O we find that spec¬ 
tral peak height of each |m) state varies depending on 
temperature, which indicates the change in thermal dis¬ 
tributions of the number states \m). At low temperatures 
(T =10 nK), the m = 2 peak is highest. On the other 
hand, at higher temperatures (T =100 nK and 50 nK), 
the m = 1 peak is dominant. To clearly discuss this 
behavior, in Fig. 01 we show the in-situ column-density 
distributions of the number states |m), which can be cal¬ 
culated by exact diagonalization for the localized Hamil¬ 
tonian in Eq. dni). At higher temperatures (T =100 nK), 
the |m = 1) state spreads widely, while larger-m states 
prefer center of the potential. Thus, Nm the total num¬ 
ber of \m) monotonically decreases with increasing m. 
This behavior is consistent with those seen in spectra in 
Fig. [3] In contrast, at low temperatures {T =10 nK), the 
Mott shell structures develop, and as a result, the num¬ 
ber state distributions for m < 2 show a dip structure 
as shown in Fig. |4l Here a three dimensional shell struc¬ 
ture is projected onto the dip structure in the column- 
density distributions. Thus, becomes nonmonotonic 
and Nm =2 becomes maximum, leading to a large inten- 
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sity of the peak of \m = 2) in spectra (see Fig. [3]). These 
features suggest that the present spectroscopy can be a 
thermometer. 

We finally compare these results with the experimen¬ 
tal observations in Ref. 0, in which spectra and also 
in-situ number-state distributions have been measured. 
We find that both spectra and distributions agree well 
with the calculations at higher temperatures than those 
at lower temperatures. It thus suggests that temperature 
will be an important parameter to discuss the quantita¬ 
tive aspect of experiments. Although nearly pure EEC 
can be created before loading atoms into a lattice, the 
loading usually induces a certain amount of heating (2^ . 
In the present calculations, a harmonic trapping poten¬ 
tial is used by focusing on the bottom of the experimental 
trapping potential 0 , whereas in Ref. 0| , by carefully 
considering the anharmonicity of trapping potential, ex¬ 
perimental spectra can be reproduced even at zero tem¬ 
perature. The calculations considering both anharmonic 
potentials and thermal fluctuations may be required for 
more quantitative agreements. 

V. SUMMARY 

In summary, we theoretically investigate the mi¬ 
crowave and the laser spectroscopy on the bosonic Hub¬ 
bard systems. We first discuss the sum rules of spectra 
up to the second order, which can be derived from the 
general principle of spectroscopy. This principle provides 
various useful information on the physical properties of 


spectra. The spectra of superfluid states with phase co¬ 
herence is broadened by the many-body effects, and its 
broadened width can be characterized by number fluctu¬ 
ations in thermal equilibrium. In contrast, the spectra of 
the number definite Mott insulating states are broadened 
by quantum fluctuations caused by tunneling effects. 

We next propose a two-mode approximation to calcu¬ 
late spectra at finite temperatures. This approximation 
assumes that spectra can be decomposed into two con¬ 
tributions originated from Bose-Einstein condensates and 
uncondensed normal states. Our method is built up by 
considering the spectral characteristics figured out from 
the sum rules, so that the multi-peak structures resulting 
from coexisting of superfluid and normal states at finite 
temperatures can be successfully dealt with. 

Finally, by combining the two-mode approximation 
with the finite temperature Gutzwiller approximation, 
we numerically calculate spectra by considering realistic 
experimental parameters of the microwave and the laser 
spectroscopy. We find that our method can reproduce the 
essential features of spectra in experiments. We also dis¬ 
cuss the lattice depth dependence and the temperature 
dependence of the microwave spectra. These results clar¬ 
ify that the present spectroscopy can be sensitive tools 
for investigating the quantum phase transitions at finite 
temperatures. 
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